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Abstract 

We address the issue of the quark-hadron duality in the spectral densities induced 
by the heavy-light quark currents Qq. In the limit uiq — > oo, m q — > we observe 
an enhancement of the physical spectral density compared to the quark one in the 
scalar and axial channels, due to the Goldstone meson contributions. This may 
imply that the scale where duality sets in is higher in these channels than in the 
vector (pseudoscalar) case. Implications for the nonleptonic decays of B mesons 
(the b — > ccs transition) are considered. We discuss in detail the decay pattern and 
obtain an independent estimate of the "wrong" sign D yield. 
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1 Introduction 



A large and rapidly growing number of applications of the heavy quark theory 
requires predictions in the Minkowski domain. The inclusive width of the heavy 
flavour hadrons is one the best known examples. All predictions of this type are 
based on the so-called QCD duality (for a recent discussion see, e.g. Ref. 

The issue of duality is as old as QCD itself. Because of its complexity it was 
virtually neglected for a long time, basically after the classical paper of Poggio et 
al. . Recently the interest to this question was revived by the necessity of having 
QCD predictions valid up to per cent accuracy in several problems of practical 
importance. The first attempts at quantifying possible deviations from duality were 
presented in Refs. [|], f|. 

It is not easy to estimate the scale of the duality violations from fundamental 
QCD. Therefore, it is of paramount importance to gain experience in various particu- 
lar problems. The present work investigates duality violations in an applied aspect. 
We will be mostly concerned with the spectral density in the two-point function 
induced by the heavy-light currents QTq, where Q and q are the light and heavy 
quarks, respectively. It turns out that in the limit vtlq — > oo, m q — > the spectral 
density in the scalar (axial) channel is very peculiar: it starts below the position of 
the ground-state resonance, and is enhanced by the Goldstone meson contribution. 
Qualitatively the emerging picture is quite different from that taking place in the 
"conventional" vector (pseudoscalar) channel, where the ground state resonance is 
physically the lowest-lying state. A key role in this phenomenon belongs to the 
spontaneous breaking of the chiral symmetry for the light quarks. 

The formulation of the problem to be discussed below is, in a sense, complemen- 
tary to Ref. 0. If in the latter work one approaches from the high-energy side (i.e. 
one considers the oscillation and asymptotic zones, in the nomenclature of Ref. [[3] , 
see Sect. 5.2), where duality is defined point-by-point; here we mostly deal with the 
resonance zone. Correspondingly, by duality we mean here that certain finite-range 
integrals over the hadronic spectral density are equal to the same integrals over the 
quark spectral density. The equality (or, rather, deviations from it) is checked by 
saturating the hadronic contribution by exclusive modes (cf. Ref. 0]). The size of 
an integration interval large enough for duality to take place characterizes the scale 
of structures and disturbances due to low-lying resonances. 

In the second part of the paper, as a practical application we consider a partic- 
ular physical process, namely b — > ccs. Due to a relatively small energy release in 
this transition it is a natural suspect || ||, [7| for duality violations, which might be 
relevant in the so-called "semileptonic branching ratio versus n c problem" ||. The 
model of Ref. estimates possible violations of duality in the b — > ccs only at the 
per cent level. If so, they can be neglected, and the duality-based predictions || [TIJ 
should be valid. The above model is admittedly crude, however, and independent 
estimates are badly needed. Since we observe a strong enhancement of the spectral 
density in the axial cs channel, associated with the soft kaon contribution, a corre- 
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sponding enhancement in the B — > D*D*K could be naturally expected. A priori, 
it could then be duality violating. 

We suggest to use factorization, and the so-called generalized small velocity 
(GSV) limit in the b —>■ ccs transition. Both approximations are described at length 
below; here we just note that using them opens all decays of the type B — > DD + 
S'-wave kaons (kaons and pions) for applications of the soft pion technique fllT 



The most economic way to implement the idea is by using the effective Lagrangian 



approach, which combines the chiral and heavy-quark symmetries fl2"|, [TBI [14 



Surprisingly, our estimates of the decay modes B — > D^D^K, where the kaon 
is treated as a soft Goldstone boson, give no indications on significant duality vio- 
lations. We hasten to add, though, that the accuracy of the above estimates is not 
high. 

The mechanism we discuss does not differentiate noticeably between the B meson 
and A& baryon decays. Thus, it adds nothing new to the problem of the A& lifetime. 
However, the l/A^-motivated factorization used in the B decays may no longer be 
applicable when a heavy baryon is produced in the final state of A& decays. This 
may provide a mechanism for the observed lifetime difference, but we will not go 
into this issue here. 

An interesting lesson follows for the lattice analysis of the heavy-light systems 
in the scalar and axial channels, where the quenched approximation is expected to 
give results different from those with light dynamical quarks. 

The paper consists of two distinct parts. In the first part we discuss general 
features of the spectral densities corresponding to the two-point functions induced 
by Qq in the limit uiq —>■ oo, m q —>■ 0. (Sect. 2). Implications for the transition 
b — ► ccs are considered in the second part (Sects. 3 and 4). Our conclusions are 
summarized in Sect. 5. 



2 Spectral densities and duality in the scalar and 
pseudoscalar Qq channels 

As a preparatory step let us consider the spectral densities for the two-point func- 
tions induced by the current Qjfj,(l — 75)9- Since our consideration at this stage 
focuses on qualitative aspects, it is convenient to work in the limit where the mass 
of the heavy quark mq — > 00 while that of the light quark m q — > 0. In the chiral 
limit all mesons (ir, K, rj) belonging to the Goldstone octet are massless. The effects 
due to the finite quark masses will be incorporated later, when we pass to b — > ccs. 
The vector and pseudoscalar Qq mesons are degenerate. The same is valid for the 
scalar and axial ones. The splitting between the first and the second multiplets is 
of order Aq CD (practically, about 500 MeV). 
The first analysis of the two-point functions 

n2? = i J Ae" 1 (o|y„(rr)yJ(o)|o> = + ... , 
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nj# = • J d 4 xe^ (o\M*)> 4(o)|o) = + ... , 

U^=i Jd A xe iqx (0\S(x),S^O)\0) , 

n( p )=2 J d 4 xe igx (0|P(a;),P t (0)|0) , (1) 

dates back to the early days of QCD [|l^]. The currents V^, A^, S and P are defined 

as 

V M = Ql^q, = Q7 M 7 5 g, S = Qq, P = Qi"f 5 q . (2) 

Although the weak currents are V — A, the scalar and pseudoscalar two-point func- 
tions appear as their longitudinal parts, 

S= — , P = d,A^ . 

m Q -m q rriQ + m q 



Already in Ref. |15[ it was noted that the quark condensate correction is en- 
hanced in the heavy-light quark case compared to, say, the classical p- meson sum 
rule |I6| , producing stronger disturbances of the spectral density at small energies 



and, thus, leading to stronger deviations from duality []. First of all, its dimension 
in the former case is 3, while in the latter case it is 6. Second, the numerical co- 



efficient is larger. Further studies [17, 18, 19] confirmed this observation in a more 



quantitative way. Actually, Ref. (T7j was the first to introduce the heavy quark 
limit. It was noted there that in this limit the spectral densities in the (transverse) 
vector and pseudoscalar channels are degenerate. The same is valid for the scalar 
and (transverse) axial channels. For this reason we will limit ourselves in this section 
only to S and P channels. 



Even more remarkable is the fact |17j that using the standard "first resonance 



plus parton-like continuum" model of the spectral density produces an abnormally 
large residue of the ground-state resonance in the scalar (axial) channels, an order of 
magnitude larger than that in the pseudoscalar (vector) channels. This observation, 
which went unnoticed, must be viewed as a hint that something unusual takes place 
in the scalar channel. Now we are ready to explain this anomaly. 

As a matter of fact, the standard model of the spectral density mentioned above 
does not work in the scalar (axial) channels because of an enhancement at low 
energies due to the contribution of the Dn intermediate state. (The mesons of the 
type D and D* built from Qq and degenerate in the limit mq — > oo, m q — > will be 
generically referred to as D's. The massless pseudoscalar mesons will be generically 
referred to as 7r's). 

In the free-quark approximation the four spectral densities (for lT y ' p ) and lT 5 ' A )) 
are all identical and equal to 

Imn=^ 2 ; (3) 

1 We stress again that we focus now on the resonance zone, where a typical scale is set by low- 
dimension condensates, unlike Ref. || where the emphasis was on the oscillation and asymptotic 
zones, and on the modelling of the impact of high-dimension condensates. 



3 



only chirally odd condensate corrections distinguish between the first and the second 
pairs of above currents. Here e is the excitation energy measured from the quark 
threshold (i.e. from the heavy quark mass mq). Logarithmic factors due to anoma- 
lous dimensions and a s corrections are neglected in Eq. @. Although they play a 
role numerically, see e.g. Ref. [0, qualitatively they are not important. Equation 
(0) will set a natural scale for the spectral density in the heavy-light channels. Its 
important feature is a strong suppression at low energies, as e 2 . 

Now, let us turn to the physical (hadron-saturated) spectral densities. In the 
pseudoscalar channel we encounter a more or less familiar picture. The physical 
spectral density starts from the ground-state pseudoscalar D, then there is a gap, 
and then multiparticle intermediate states (continuum) add up into a curve that is 
expected to be relatively close to the quark one. The continuum actually is a sum 
over broad resonances, the radial excitations of D. 

The first resonance in the pseudoscalar (vector) channel is situated at e = A 
where A = Mo — m c \ the continuum threshold is at e = e c ~ 2A, see Fig. 1. The 
onset of "continuum" approximately coincides with the position of the first radial 
excitation in the given channel. The parameter A above is a basic parameter of the 



heavy quark theory @. Numerically A ~ 500 MeV [|T], ||, |3[. (We recall that 
although we keep using the name U D mesons" so far we work in the limit mq — ► oo.) 

The scalar (axial) resonances we are interested in are P-wave in the language 
of the naive quark model and l/2 + in the language of HQET. Due to the chiral 
symmetry breaking they lie higher.0 (For a concise review of the higher Qq states 
see Ref. pi| .) The heavy quark symmetry implies that the scalar Dq and the axial 
D\ are degenerate; one expects both Dq and D\ at e ~ 2A. What is remarkable is 
that the scalar (transverse axial) spectral density is non-vanishing and large below 
Dq (Di), in the interval from A to ~ 2A, since it receives an enhanced contribution 
from the states Dtt (in the scalar channel) and D*tc (in the axial channel). The 
pions are strongly coupled to the Qq current, in the S wave. 

Consider for definiteness U^ s \ In the low-energy limit the pion is soft and can 
be reduced by using the soft-pion technique. Then, 

(0\Qq\Dn)^Q = l( |Qi 75 g|L>) = -&m Q . (4) 

Jtt Jit 

Here fq is defined as 

{®\Qlnl5<l\D) = ifqPf,, 

and we do not differentiate between Mo and mq in the limit mq — > oo; k is the 
pion momentum. The spectral density takes the form 

Imn< s ) = l -MQq\D*)\ 2 x phase space = (ty - ^ ) (e - A) , (5) 

2 The scalar resonance will be referred to as Dq. The axial l/2 + resonances are sometimes 
denoted as D[, to distinguish them from the axial 3/2 + resonances denoted as D\. Since in the 
present paper we will discuss only the axial l/2 + resonances, the prime will be omitted. The 
primes are reserved for radial excitations. 
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Figure 1: A sketch of the spectral density versus e/A in the pseudoscalar channel. 
To set the scale we show also the free quark spectral density represented by the 
parabola. 
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where Nf = 3 is the number of massless flavours (Nj — 1 is the number of Goldstone 
mesons) and e — A is the energy measured from the position of the ground state D. 

Comparing Eqs. (|3|) and @ reveals the remarkable enhancement of the spectral 
density at low energies that was mentioned above. First, the dependence on e is 
parametrically different. Due to the S'-wave nature of the matrix element (0\S\Dtt) 
the Dir contribution to the spectral density vanishes as 0(e) while the free quark 
expression (^) is 0(e 2 ). Second, as is well known, the static coupling fq is rather 
large numerically |T^, p5|| . The combination f^rriQ stays constant in the limit niq — > 
oo [IT]], modulo the hybrid logarithms (EBfl, and this constant is approximately 0.2 
to 0.25 GeV 3 g, g5[. Treating Eq. (|) literally we get that the Dn contribution to 
the spectral density at e = 2 A exceeds the free quark expression by a factor of ~ 1.5 
(Fig. 2). The integral over the Dir spectral density up to e = 2 A is approximately 
1.2 times the integral over the free quark spectral density from to e = 2 A. 

Of course, at a certain energy, the soft-pion result (f|) becomes invalid - the 
momentum-dependent part of the interaction will cut off the amplitude. The cut off 
presumably occurs near the position of the resonance D$. The excess of the spectral 
density below Dq and at -Dq has to be compensated by an extended gap immediately 
after Dq. This gap stretches, presumably, up to e ~ 3A. In general, the situation is 
very similar to what occurs with the two-pion contribution in the spectral density 
induced by the gluonic current a s G 2 ||27|1 , where a strong low-energy enhancement 



is supplemented by an extended gap, and the scale of the duality violation is quite 
large. 

Note that the contribution we have calculated has proper dependence on 
parameters N c and Nf. In the limit N c — > oo it decouples compared to Eq. (§), 
as it should. There is no decoupling, however, if Nf/N c stays fixed. The pattern 
we observe presents another example, similar to Ref. pTj, that in the low-energy 
domain the 1/N C counting rules should be taken with caution. There exist certain 
mechanisms that can totally upset the 1/N C estimates at N c = 3 (although at 
academically large N c the 1/N C counting will work, of course). 

Thus, in the pseudoscalar channel we expect the duality interval to be ~ 2A, 
while in the scalar channel its size is expected to be larger, ~ 3A. By the duality 
interval we mean the following: the minimal interval of energy (0, £q) needed to 
make the smeared resonance contribution (approximately) equal to the quark-gluon 
one. In other words, in the smearing integrals over the spectral densities from to 
Eo the upper limit Eo should be chosen at ~ 2A in the pseudoscalar channel and 
~ 3A in the scalar one. 

We pause here to make a few comments. First, if we consider the sum of the scalar 
plus pseudoscalar heavy-light channels (or vector plus axial) the quark condensate 
correction (qq) cancels, leaving us with no hint that the duality violation scale is 
larger in this case than in the classical p meson sum rule, at least at the level 
of analyzing the first terms in the operator product expansion. This fact is rather 
alarming, since in the weak inclusive decays we deal with the V — A currents, and the 
above cancellation is quite typical. The model of Ref. |3f gives no clues in this case 
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either, since it admittedly omits all effects specifically related to the spontaneous 
breaking of the chiral symmetry. At an observational level we detect here a possible 
signal that the scale of the duality violations for V — A currents may be larger than 
can be inferred from the analysis of the lowest-dimension condensates. Does this 
mean that some other condensates, of higher dimension, must be important? And 
how can they be identified? These questions still remain open. It is clear that the 
sum rules for the heavy-light currents have to be reanalyzed with emphasis on this 
aspect. In particular, an updated analysis of the scalar (axial) channel is more than 
in order. One has to replace the standard "lowest resonance plus continuum" pattern 
by the more complicated picture described above. The impact seems to be obvious. 
The residue of the Dq (D{) state will go down with respect to the prediction of Ref. 
T7| . A part of the spectral density will be pumped out from the resonance into the 



non-resonant Dtt background. 

The second remark refers to the lattice calculations of the heavy-light systems 
(spectra and coupling constants). Pseudoscalar and scalar channels should be dras- 
tically different with respect to inclusion of the light dynamical quarks. If the 
pseudoscalar channel seems to be relatively stable, unquenching quarks and making 
light quarks really light must produce a dramatic effect in the scalar one. Of course, 
this is not the first time the dynamic light quarks lower the threshold. For instance, 
in the p meson case there is the two-pion cut as well. This state, however, is only 
relatively weakly coupled to the current q^f^q, since the pions are in the P wave; 
therefore, speaking in practical terms, this contribution is rather unimportant (al- 
though at asymptotically large separations it will dominate anyway). This is not the 
case for the current Qq, where the intermediate state is not only lower in mass 
than Dq (Di), but also strongly coupled to the current. Therefore, its impact in the 
two-point function should be essential and the correlator must drastically change 
once the light quarks are unquenched and made light. Studying this problem on the 
lattice seems to be a nice testing ground for various approximations routinely made 
within this approach. 



3 Nonleptonic B decays (6 — > ccs): outlining the 
problem 

Having established the enhancement of the Dn intermediate state in IulIT^ 5 ) (or 
D*n in near the threshold, it is natural to turn to nonleptonic decays of 

B mesons, the b —>■ ccs transition. Indeed, in this case the cT^s current produces 
charmed/strange hadronic states with the quantum numbers 1~, + (the vector 
current) and 1 + , 0~ (the axial current). The 1 + and + channels were shown above 
to be responsible for an enhanced production of Goldstone mesons in the S wave. 
Of course, in the actual B decays the situation is not quite the same as in the limit 
ttiq — > oo , m q — > considered in Sect. 2. Suffice it to mention that the actual 
values of the c and s quark masses are such that the axial ground-state meson D sl) 
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shown in Fig. 2 at the end of the shoulder (at e ~ 2A), turns out to be almost 
exactly at its beginning, presumably barely above the threshold of D*K. The l/2 + 
charmed strange mesons have not been detected experimentally so far. Moreover, 
the actual value of j\ is also noticeably lower than its static value. Also, the energy 
release in the light particles in the case at hand is < 1.5 GeV. It is quite clear that 
the kaon mass is not negligible in the phase space. For this reason in the nonleptonic 
B decays a dedicated analysis is needed. 

First, however, a few remarks regarding the general situation with the b — ► 
ccs transition are in order. This transition came under renewed scrutiny recently 
[7[ on purely phenomenological grounds, in an attempt to find a solution of the 
"branching ratio versus n c problem" . In Ref. it was assumed that the theoretical 
understanding of b — > cud is solid. Then, from the measured semileptonic branching 
ratio of approximately 10.5%, the branching ratio of b — > ccs was predicted to be 
close to 30%, with the corresponding charm yield n c ~ 1.3. Moreover, using the 
most naive duality estimates in conjunction with the parton model for b —>■ ccs, it 
was suggested @ that approximately 1/2 of b — > ccs hadronizes in the channels of 
the type B -> D^D^KX, and only the remaining 1/2 goes to D^D^X, the 
channels on which the attention was focused previously. If so, about 15% of the B 
decays have to end up with the "wrong" sign D's and K's in the final state. 

This prediction was confirmed by recent results from the CLEO, ALEPH and 
DELPHI experiments |]28]| , reporting the yield of the "wrong" sign D's at the 10% 
level. The corresponding value of n c is close to 1.24 f28|j . Additionally, ALEPH has 



recently reported a value for n c in Z — > bb [291 of 1.23 ± 0.07. 

Although the situation can be looked at quite optimistically, serious questions 
are still to be answered. The first question is, of course, purely experimental. As 
was emphasized in Ref. ||28||, the charm yield n c , as computed in the usual way from 



the measurements at T(4S), remains unchanged: n c = 1.10 ±0.06 |30|]. The contra- 
diction is obvious, suggesting that the experimental situation is still not settled. 

We clearly cannot comment on this aspect, and quickly pass to what is known 
theoretically. Since the leading nonperturbative corrections are expected to play no 
significant role in the issue ||, the focus of theoretical analysis is shifted towards 
perturbative calculations in b — > ccs. The first dedicated calculation of the gluon 
corrections was carried out in Ref. |3l . The most advanced analysis existing today 



is presented in Refs. || [T(J (see also references therein). To quote representative 



values of the predicted parameters, let us assume that a s (Mz) = 0.11 and fi = m&/2, 
in MS; then | n c = 1.28 (the corresponding BR s i(i?) is slightly lower than 11%). 
Within a somewhat different procedure of treating the ratio of the quark masses 
rric/rnb (but the same values of the parameters as above), the theoretical numbers 
become P2J n c = 1.23 and BR s i(5) = 11.5%. Thus, it is fair to say that at the 
current level of understanding the theoretical prediction for n c is close to 1.25. 
Without the 0(a s ) correction the charm yield is 1.15. Thus, the inclusion of the 
0(a s ) gluon correction enhances BR(6 — * ccs) by a factor ~ 1.6. 

In the approximation of Ref. || factorization is explicit. This statement will be 
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explained in more detail below (Sect. 3), where we will introduce the factorization 
hypothesis, one of the key elements of our consideration. The impatient reader 
may turn to Ref. |I3| for very clear explanations. Here we just note that, in any 
perturbative calculation respecting factorization, the transitions corresponding to 
the vector cj^s and the axial 07^753 currents necessarily have equal probabilities, 
provided very small effects cx m 2 s are neglected. If so, the result || implies that the 
branching ratios of b — > ccs, with cs from the vector and axial currents, respectively, 
are 12 to 13% eachfj 

It is clear that the theoretical predictions || discussed above are quite com- 
patible with the newest experimental trend. The basic assumption underlying the 
theoretical analysis is the validity of the quark-hadron duality. Strong violations of 
duality in b — > ccs would destroy the predictive power of the existing theory. Al- 
though at the moment no sources for such large violations were identified ||], in the 
absence of a complete theory it is obviously desirable to have as many independent 
confirmations as possible. We will analyze below the transition b — > ccs from this 
point of view. 

The transition b — > ccs is singled out in this aspect for the following reason. The 
relative smallness of the energy release, which alerts one regarding duality violations, 
can be turned into an advantage. Indeed, in this case the number of the hadronic 
channels saturating the physical decays is not large, and one can try to estimate 
these channels one by one to see whether they add up to the quark-gluon result. As 
a by-product one could hope to get a more direct estimate for the B -> D^D^KX 
rate. The goal here is to check whether the 10% yield reported experimentally is well 
understood theoretically, invoking as few unsubstantiated assumptions as possible. 

4 Nonleptonic B decays (b —> ccs): analysis of ex- 
clusive modes 

The relevant part of the weak Lagrangian contains two structures, with the "direct" 
and "twisted" colour flow 



L = ^v d ;\\, 



V2 



ai(6r M c)(cr M s) + a 2 {bT (t s){cr (t c) . (6) 



We will disregard the second term, with the twisted color structure, for the following 
reason. First, the value of I02/ oti| is rather small numerically (see e.g. the review 
paper P0fl), approximately 0.1. Therefore, the square of the second term contributes 
at the level of 0.01. The interference with the first term is suppressed by N c and 
is thus expected to show up at the level of corrections ~ 0.2/3 in the probability 

3 In Ref. |3J| some of the 0{a 2 s ) corrections were estimated, and at this level nonfactorizablc 
terms appear. For the purpose of our analysis one can safely use factorization, at least as a starting 
point. 
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of the decay modes we are interested in. The estimates to be presented below have 
intrinsic theoretical uncertainties of this order of magnitude or larger. 

Furthermore, in considering the first term it is reasonable to accept, at least at 
this stage of the analysis, the factorization approximation. By this we mean that in 
treating the hadronic matrix elements, the bT^c bracket of the effective Lagrangian, 
will be factored out from the cT^s bracket. The corresponding hadronic subsystems 
are assumed not to communicate with each other through the soft gluon exchanges, 
although inside the subsystems all these exchanges are taken into account in full. 
Note, that we also automatically include all hard gluons (with off-shellness larger 
than rrib), through the factor a%. The standard theoretical justification for the 
factorization hypothesis is the 1/N C counting. All non-factorizable contributions are 
suppressed by 1/N C . Note that the modes with the hidden charm (e.g. J/ipK) will 
not concern us here. 

Certainly, we are well aware that the non-factorizable contributions must be 
present (see e.g. ||34| ) ; their effect is noticeable in the fine structure of the nonleptonic 
decays and in some special modes, but otherwise it is quite modest. For instance, in 
B — ► DD S the non-factorizable contributions were estimated to be less then 10% of 



the factorizable part p5 ]. We will ignore the non-factorizable terms in the present 
work. 

Another theoretical tool which will help us is the generalized small velocity limit, 
i.e. the assumption that two charmed mesons we deal with in the b — > ccs transition, 
in the final state, are slow. Kinematically this means that 

M B - 2M D < M D . (7) 

This is not the first time the GSV limit is exploited in the context of B decays, see 



e.g. [|3j|, [I(J. Using the GSV limit, in combination with factorization, will allow us to 
disregard excitations in the bT^c bracket, and exploit the well developed formalism 
of the soft Goldstone mesons for the modes of the type DDK. 

Although in purely theoretical aspect the GSV limit is an excellent tool, in the 
actual B decays Eq. (0) is valid only marginally. We are neither too close to 
the GSV point nor too far from it (cf. e.g. p3||). Under the circumstances, in 



kinematical factors we will keep the terms containing the D meson velocities, while 
omitting v 2 terms when they are additionally suppressed. In the future, with more 
phenomenological input, this approximation can be improved, including all terms 
quadratic in the (spatial) velocities. 

To warm up let us consider two-particle decays. This exercise is not new (see 
Ref. 0j), and we repeat it merely to introduce our notation and explain the choice 
of numerical values of the relevant parameters. 

Thus, factorization and the GSV limit are starting elements of our analysis. If so, 
then the bracket bT^c is responsible for the B — > D (B — > D*) transition. Production 
of excitations by this bracket - either resonances (say, radial excitations of D and 
D*), or nonresonant states of the type Dn - is suppressed by the velocity squared 
of the charmed meson. The Dtt production by bT^c can be studied within the chiral 



10 



perturbation theory |38[ . One has to consider pole graphs which, in addition to the 
velocity suppression, are proportional to the D*Dtc coupling constant g. The latter 
was calculated within the QCD sum rules [39|, and turns out to be rather small, 
g ~ 0.3. Therefore, here and below we will consistently disregard all contributions 
to the decay rate that are proportional to v 2 and g 2 . 

Limiting ourselves to the D and D* states in the bT^c bracket gives us, by itself, 
a predictive power since the form factors of the B — > transitions at zero recoil 



are normalized to unity |^0|, ^T[, and near zero recoil are well approximated by the 
first derivative of the Isgur-Wise function £ [fl^| . 

The second bracket, cT^s, creates sometimes D s or D* states, sometimes radial 
excitations of D s and D*, and sometimes nonresonant D^*'K pairs. The axial cur- 
rent, in its transverse part, can produce D sl and excitations, while the longitudinal 
part of the vector current can produce D* and excitations. Let us first concentrate 
on D s and D*. 

(i) B —* D^*'Dy\ "Wrong" spin correlations 

The amplitudes for two-body transitions are given by 



A(B DD S ) = ^V cb V* s ai ( 2^M B M D f D M D 



M B - M D 1 + w' 



M 



D s 



(8) 



A{B - D*D* S ) = ^V cb V: sai ^M B M D ,f D M D ^ x 



V2 

1 ?' 
--{ev){ev) + -e^p\e^e v v p v x 



(9) 



A(B^DD:) = ^V cb V: sai x 



2JM B M D f Dt M Ds 



M B + M D 



2M, 



D 



e"v) £{vv') 



(10) 



A(B^D*D s ) = ^=V cb V; sai x 
M B + M D , 



2y/M B M D .f D .M Ds 



2M 



D, 



Here we omitted inessential overall phase factors in front of some amplitudes. The 
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constants fp s are defined as , 

(D a \c W s\0) = -if Ds M DeV ; , (D* s \c lfl s\0) = f Ds M Dt e; . 

The pseudoscalar and vector constants are identical in the limit m c — > oo, but may 
differ a little, due to preasymptotic terms, for the actual c quarks. For simplicity 
the difference between the pseudoscalar and vector constants are neglected within 
any hyperfine multiplet. Anyway, we do not know them to that accuracy. Moreover, 
v is the four-velocity of the decaying B meson, v' is the four-velocity of the 
meson and v" is the four- velocity of D^*\ 

, Ml + Ml- M* „ Ml + Ml s -Ml , „ Ml -Ml- M* a 

vv = , V v = 2 , V V = 1 , 

2M B M D 2M b M Db 2M d M Ds 

and likewise for other decays, e' and e" are the polarization vectors of D* and D*, 
respectively. 

In the limit of slow D's, when vv', vv", v'v" —>■ 1, these expressions essentially 
simplify. Thus, in Eq. (§) the square bracket becomes unity, and in Eqs. (|), ( |T0| ) 
and ([□]) the square brackets tend to zero - as yjvv" — 1 in Eq. (jT^) and as yjvv' — 1 
in Eq. (pTTJ) . Thus, in this limit we have a rigid spin correlation. If D is produced by 
the be bracket, the cs bracket will produce a pseudoscalar D s ; if D* is produced by 
the be bracket we get, in association, a vector D*. This is because in this limit the 
B — > D transition is caused by the zeroth (time) component of the current, whereas 
B — > D* is due to the spatial component. The "wrong" -spin transitions, Eqs. fllPP 
and fllTD , switch off in this limit. Table 1 gives the values of the parameters y = vv', 
v'v" and vv" for various transitions in the ground state and excited resonances. In 
the GSV limit all these kinematic parameters reduce to unity. 

This table also gives the values of "block factors" . They are defined as 

block factor = — 7 — ^ = ^ = 

\\^V cb V* ai (2y/M^MEf Da M Da 

where the summation over the polarizations of Df* is implied, where applicable; 
for the excited D s states considered below Md s and fo s in the denominator will 
be understood as the mass and the coupling of the respective pseudoscalar. We 
included here the IW function see below. In the first decay the block factor 

is just the square bracket in Eq. @. 

We see that in the B — > DD S transition the block factor which is equal to unity 
in the GSV limit, is actually close to 1.5, while in the "wrong" spin correlation 
transition B — > DD* the block factor, which vanishes in the GSV limit is actually 
close to 1.25. This is a manifestation of the fact that in the two-body decays into 
the lowest-lying D's the GSV approximation is not too good, as was expected, of 
course. 

Table 1 shows that the GSV approximation becomes much better for excited .D's 
where the wrong spin transitions are indeed suppressed. The GSV approximation, 
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Table 1: 



Decay 


V 


v'v" 


vv" 


block factor 


B _> DD S 


1.39 


2.79 


1.36 


1.50 


,S 


1.36 


2.53 


1.29 


1.25 


B -> B*B S 


1.32 


2.52 


1.33 


.93 


B -> B*B* 


1.29 


2.28 


1.27 


^■0.94 


b -> dd; 


1.25 


1.82 


1.13 


1.22 


B -> BB*' 


1.21 


1.64 


1.10 


0.79 


B -> B*B^ 


1.19 


1.64 


1.11 


0.58 


B -> D*D*' 


1.15 


1.47 


1.08 


1.04 


B -> BB" 


1.10 


1.28 


1.04 


1.07 


B -> BB*' 


1.05 


1.14 


1.02 


0.37 


B -> B*B^' 


1.05 


1.14 


1.02 


0.28 


B -> D*D*' 


1.006 


1.02 


1.002 


v/3- 1.09 



quite naturally, significantly improves in the three-particle decays as well, where a 
part of the overall energy release goes to create, additionally, the kaon mass, and 
the remainder is shared by three particles, not two. 

If instead of D s and D* we have their radial excitations, Eqs. ©-(01) are 
modified in a minimal way. Apart from the masses and the kinematical factors, one 
must replace fu a by the corresponding coupling. 

For numerical estimates of the decay rates we need to fix various parameters. 
First, we take f Ds w 200 MeV fE|. Now, the Isgur-Wise function £ must be 
evaluated at the proper recoil values, which are given in Table 1. We use the linear 
approximation for the Isgur-Wise function 

£(W) = l-p 2 (W-l), (12) 

where p 2 is the slope parameter. Its numerical value is more or less known, and we 
put 

p 2 = 0.7. 



This value is compatible with the experimental data [0 and with the QCD sum 
rule calculations |4"5fl . Finally, w 1.1. 
The relevant decay rates are 

r(B ^ B ,. ro ("i±^)V, (13) 



T(B -> DDI) 



{M B + M D ) 
AM 2 D 



(14) 
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Table 2: 



Decay 


Branching , % 


B -> DD S 


1.1 


B -> DD* 


0.87 


B -> £>*£> s 


0.45 


£ -> 


1.5 


Total Br, % 


4.0 



F(B -> 



(Mg + M D ,) 2 
r ° 4MB (( } 



(15) 



Here 



3r£ 



1 + 



+ 



A 2 



W) 2 - 1 
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+ 



(w')(w") 



e 2 - 
J (16) 



G 2 



i o = -£\V A V„\*a\ f Ds M D M 2 Ds ^-, 

and r *, T' Q , Tq are the same with the D (D s ) replaced by D* (D*) for asterisk, 
and p is the c.o.m. momentum of the produced mesons. Numerically the relevant 
branching ratios are given in Table 2. 

So far we basically repeated the calculations of Mannel et al. ||37|| ; our estimate 
of the Isgur-Wise function is different, and our value of fo s is essentially lower than 
that accepted in Ref. |J7|. We note that the total sum of the above four modes 
amounts to ~ 4% in the branching ratio. None of these modes produces "wrong" 
sign D's. The above numerical results are in agreement with the experimental data 



T(B° — > DD S ) = 0.7 ± 0.4% , T(B° — > D*D* S ) = 1.9 ± 1.2% , (17) 

Br(B° -> D*D S ) = 1.2 ± 0.6% , Br(5° -> £>D*) = 2 ± 1.5% , (18) 

Br(fi+ -> D5 S ) = 1.7 ±0.6% , Br(5+ -»■ D*D* S ) = 2.3 ± 1.4% , (19) 

Br(£° -> D*D S ) = 1 ± 0.7% , Br( J B° -> DL>*) = 1.2 ± 1%) . (20) 

Encouraged by this success we can proceed to other exclusive modes which have 
not been discussed in the literature so far. 

(ii) Modes with the radial excitations of Dy> 

We got used to the fact that the radial excitations are coupled to corresponding 
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currents significantly weaker than the ground-state mesons with the given quan- 
tum numbers. This is, for instance, the case with the J/i/j and T mesons whose 
coupling to the current Q^^Q decreases with the excitation number. In the decays 
B -> £>(*)£>(*)' (where the prime denotes the first radial excitation and the dou- 
ble prime will denote the second excitation) we expect the opposite trend. This 
is a specific feature of the Qq current. The corresponding parton spectral den- 
sity grows quadratically with the energy e, see Sect. 2. If duality takes place, 
the residues should be roughly proportional to the corresponding areas under the 
parabola in Fig. 1. To get an idea of the excited state contribution we will assume 
that M D / « 2.6 GeV, M D u « 3.1 GeV, M Df « 2.75 GeV and M D «, « 3.25 GeV, i.e. 
are equidistant. Imposing the duality condition in the form 



/ dMImn pcrt (M) = / dMlmU k (M), TmH k (M) oc M k f%6(M - M k 



we then obtain 



flap's _ 27-8 _ fp^Mp 'j . 64-27 

P D M Ds 8 ~ ' ' f D M D 



2 T\/r_ o ^2 U o V J 



Note that we use here an "extended" definition of the excited resonances, which 
need not coincide with the standard Breit-Wigner peaks, with the backgrounds sub- 
tracted. Rather, we include backgrounds, collecting, say, in D' s everything with the 
invariant mass lower than w 2.3 GeV (except the ground state D s , of course). This 
is a natural operational definition from the point of view of the QCD practitioner. 

The estimates in Eq. fl2~T|) are crude and intended only for the purpose of orienta- 
tion. However, it is a model for the spectral density that exactly respects duality to 
the given approximation, and, thus, allows simultaneous estimates of its violations 
in the actual decay width, where only a limited energy release is available. The 
contribution of the last open multiplet is actually a reasonable estimate of the scale 
of duality violations that might occur. On the practical side, to allow a measure of 
theoretical uncertainty in the excited state residues, we will allow the ratios to float, 

(f 2 D ,M D ,)/(f D M Ds ) = 2.5 x, (f 2 D rjM D ,,)/(f D M Ds ) = 4.5* (22) 

where x will be varied between 0.6 and 1. The values of x < 1 reflect a relative 
enhancement of the lower end of the spectral density by the perturbative and, to 
some extent, condensate effects. Also the higher end of the spectral density is 
somewhat suppressed since at these energies relativistic effects become important, 
and temper the e 2 growth of the spectral density characteristic to the non-relativistic 
limit. As we will see later on, the data are more consistent with x ~ 0.6. 

In Table 3 the branching ratios for the two-particle decays into excited states are 
quoted for three different values of x . Altogether we get that 

Bi(B -> D { * ] + excited D^) ~ 7 to 12% . 
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Table 3: 



Decay 


Br % fx = 6) 


Br % (x = 8) 


Br % fx = 1) 


B -> DD' S 


1.1 


1.5 


1.9 


B -> DD*' 


0.5 


0.65 


0.8 


b -> b*b; 


0.25 


0.35 


0.4 


B -> D*D*J 


2.4 


3.2 


4.0 


Total £>W' Br, % 


4.3 


5.7 


7.1 


b -> bb; 7 


1.2 


1.6 


2.0 


B -> BB*" 


0.1 


0.15 


0.2 


B -> B*B^' 


0.06 


0.08 


0.1 


B -> D*D*" 


1.1 


1.5 


1.9 


Total B w " Br, % 


2.5 


3.3 


4.1 


Total £>W excit Br, % 


6.8 


9.0 


11.3 



Let us discuss the decay pattern of these excitations. This issue is interesting not 
only in connection with the "wrong" sign B's in the B decays, but by itself as well. 

D' s can decay neither into DK nor into D*K. The second mode is presumably 
below the threshold. Even if it is slightly above the threshold, the small energy 
release, in conjunction with the P wave nature of the decay, strongly suppresses it. 
Thus, the dominant mode is expected to be D s tttt. Therefore, the D' s production 
does not generate the "wrong" sign B's. 

On the other hand, D" will presumably decay predominantly into D*K. One 
third of the above goes into D*i]. Some DK* are also possible. The latter mode 
also leads to a "wrong" sign B. Its presence may somewhat distort, however, our 
estimate of D*r]. A competition from the D s tttt mode is generally possible, but 
unlikely to be essential: the D" — > D s -n-n is a three-body decay, with the additional 
suppression due to the pion momenta. (It can proceed via an S-wave nn resonance; 
however, the a meson is too broad and leads only to a moderate enhancement, 
whereas higher resonances do not have much of a phase space.) 

As for D*' and B*", their dominant decay mode is likely to be D^K and D^rj. 
Thus, in the three latter cases we end up with the "wrong" sign B's. Due to S , B(3)a 
symmetry the yield of r) is three times smaller than that of K's. The total branching 
of the corresponding modes is 5 to 9%. The DD s rj yield from this mechanism is 1.2 
to 2.3%. 

(Hi) Nonresonant D^D^K modes 

So far we discussed the 1~ and 0~ channels of the cs current. Now let us proceed 
to the positive-parity channels. 

As is clear from Sect. 2 the nonresonant S wave kaons are produced by 07^755 
in association with D^*K The longitudinal part of the current produces DK, while 
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the transverse component produces D*K. We will assume that the kaons can be 
treated as soft Goldstone mesons. Needless to say, factorization is implied too. 
Then, we have to write out the currents bT^c and cT^s in the chiral/heavy 



quark theory. The second current is given in Ref. [12|. The amplitude of the kaon 



"bremsstrahlung" in the transition B — > D^D^K consists of two parts - contact 



and pole (see e.g. graphs of Figs. 1 and 2, respectively, in Ref. Ul2| ; see also flSql ). It 
is not difficult to check that all pole graphs vanish in the GSV limit, when Z)'s spatial 
velocities are set equal to zero. Additionally they are numerically suppressed by g, 
the D*Dtt coupling constant, which, according to Ref. [3{|, is about 0.3. Therefore, 



at the current level of accuracy, we prefer to omit the pole graphs altogether. In the 
future, as confidence in the numerical value of g grows and estimates of the excited 
meson couplings acquire better accuracy, it will be necessary to return to the issue 
and include the pole graphs in the theoretical predictions. 

With the pole graphs discarded, calculating the decay rates of the above tran- 
sitions becomes trivial. Indeed, the ratio of the amplitude D^'D^K in the limit 
Pk — ► to the corresponding two-body amplitude is given by i/ f n - (Note 

that we will consistently neglect all SU(3)a breaking effects everywhere, except the 
phase spaces. Accounting for SU(3)a breaking effects in the amplitudes in the first 
order in m s is possible; we will defer this exercise as well.) Since the kaons are 
produced in the S wave, and the contact amplitudes depend neither on the kaon 
energy nor on the angles, the ratio of probabilities is given merely by the ratio of 
the three-body to two-body phase spaces V^, 



T(B -> D^D^K) 1 V 3 



(23) 



The three-particle phase space is conveniently written out e.g. in Ref. |47j . 

Note that the above estimate is valid only for the "right" spin correlation modes, 
DDK and D*D*K. The "wrong" spin correlation modes vanish in the GSV limit, for 
the same reason as was explained above, and cannot be computed without inclusion 
of the pole graphs. Their role, however, is significantly reduced in the three-particle 
decays compared to the two-particle case, since in the three-particle decays we are 
closer to the GSV limit. 

It is convenient to present the ratio V3/V2 as follows: 

V 2 32tt 2 I 2 ' K 1 

where the first factor corresponds to the ratio of the phase spaces with all massless 
final particles. The dimensionless factors I3 and I 2 reflect the finite masses: 

V 3 = Ml/ 3 /(2567r 3 ) , V 2 = ^p/M B = A . 

47T 87T 

Numerically J 3 are summarized in Table 4, while the branching ratios for three- 
particle yield for different nonresonant decay modes are given in Table 5. 
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Table 4: 



Decay 


h 


B -> DDK (rj) 


0.073 (0.068) 


B -> D*D*K (rj) 


0.039 (0.035) 


B -> DD'K (77) 


0.014 (0.002) 


B -> D*D*'K (rj) 


0.002 (0.001) 



Table 5: 



Decay 


Br , % (x = 1) 




-> DDK 


0.5 


B - 


-> D*D*K 


1 




-> DD'K 


0.3 




-> D*D*'K 


0.15 



Using the numerical values for J3 one can easily estimate the corresponding 
branching ratios: 

r(j -> DDK) ^ T(B -> D*D*K) ^ 

t(dd s ) ~ u ' 5 ' r(D*i);) ~ ' 1 j 

and 

r(S -> DD'K) T(B -> D*D*'K) . . 

r(^) ^ 02 • r (J . 0r ) 53 °' 04 • (26) 

Here we accounted for the fact that one can have two ^'s, of different charge, in 
each process at hand, and incorporated the block factors for the two-body modes. 
The relevant branching ratios, computed directly, are collected in Table 5. 

Altogether we see that the total branching for four possible nonresonant channels 
is ~ 2%. In the same approximation it is easy to estimate the yield of nonresonant 
D^ri instead of In the SU(3) fl limit it is 1/3 of the kaon yield. Taking 

account of the differences in the phase spaces we get Br(£> — > DD s i]) ~ 0.3 to 0.5%. 

Altogether the nonresonant channels give ~ 2 — 3% of the total branching ratio. 

(iv) Other modes with nonresonant Goldstone mesons 

Invoking the chiral/heavy quark technique and factorization in the same vein as 
above, it is possible to calculate amplitudes with two soft Goldstone mesons (non- 
resonant), e.g. DDKtt or DD s 7t-k. For instance, the D*Ktt state will be produced 
by the vector cry^s current, and all relevant matrix elements are already known in 
the literature |y]]. A number of pole graphs are to be taken into account. We did 
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not attempt this calculation, because the four-particle phase space in the processes 
at hand is prohibitively small. The decays with two nonresonant Goldstone mesons 
are rare. 

(v) Resonances with the quantum numbers + and 1 + in the cs channel 

Since the non-resonant contributions in the above channels turn out to be rather 
small, well below the 12% level following from the quark-gluon calculation (plus 
duality), it is natural to conclude that D s \ and D' sl (and D* and D*' Q ) play a 
significant role. The mass of D s i is expected to be ~ 2.5 GeV, and it is natural to 
expect that the mass of D' sl is close to 3 GeV. The masses of D* Q and D*' are close 
to those above. 

The estimates of the corresponding yield can be approached in a way similar 
to what has been done for the vector (pseudoscalar) channel. (The effect of the s- 
quark mass may be harder to take into account properly, since, as we have seen, the 
significant chiral kaon contribution gets suppressed, which may affect the residues 
of the resonances.) We use the following simplified estimate. The rate of D s i is 
expected to be approximately the same as that of D*' . The first excitation of D s \ 
is produced roughly with the same rate as the second excitation of D* s . The same 
is valid with respect to D* and D*' relative to the excitations of D s . 

One may expect that all these resonances decay in the D^K mode; in some 
cases K* may appear, instead of K. This will still lead to the "wrong" sign D. 
Then one gets ~ 7 to 11% depending on the value of x. 

(vi) Excited D states produced by bT^c bracket 

Above we considered only decays of B into and some radially excited D s 
states. There may also be decays into the excited states of D. The total yield for 
such decays with the radial excitations is quite small. There exist two possibilities: 
one can produce from the bT^c current either radial excitations of D^*\ or P-wave 
mesons. 

In the first case the corresponding partial widths are proportional to the square of 
the Isgur-Wise function of the transition between B and the radial excitations of D. 
However, the corresponding IW function vanishes as vv' — l. Thus the corresponding 
partial widths are suppressed by a factor (vv' — l) 2 ~ v' 4 relative to the transitions 
into £)". A direct calculation shows that if one takes the masses of the excited D 
states to lie uniformly below the masses of the excited D s states by 150 MeV, the 
corresponding recoil values of vv' — 1 spread from to 0.2. Thus the corresponding 
partial widths are suppressed by a factor of order 0.04 relative to the decays into the 
non-excited states. Our estimate would show that the corresponding contribution is 



less than 1%; although a larger effect is expected from the power corrections |58, (48 
up to 10 to 20% of the elastic transition, they still seem to be small. 

Higher rates can be due to the P-wave mesons from the bT^c current; in this case 
the rates are only linear in vv' — 1. An upper bound (and a reasonable estimate) of 
them can be obtained by merely calculating the sum over all decay modes with D 
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and D*, with a "faked" p = 0.25 (or with the excitation-free IW function £o{vv') = 
((1 + vv')/2vv') 1 / 2 ), through the sum rules for the transition amplitudes [fi8[| . We 
then compare the result obtained to that with the "actual" £ of Eq. (P~2j): 

Bi(B - D cxcit +X 5s ) « Br(S^ J DW+X es )|^ -Br(S^L>W+X gs )|^ =1 _ p2{w ,_ 1) . 

Since for higher excitations of Z}(*) the phase space is limited, this relation only 
overestimates somewhat the associated production of excited D^*\ In this way we 
arrive at the conclusion that the contribution of the excited states from the 
bT^c bracket is ^ 3%. 

5 Conclusions 

The existence of the Goldstone mesons results in a very peculiar pattern of the 
heavy-light spectral densities, which turn out to be drastically different, say, in the 
vector and axial channels. By treating kaons as soft, we are able to calculate the 
nonresonant D^D^K (D^D^tj) yield in B decays. 

Altogether, summing all two-body and three-body decays modes, we get 23 to 
31% for the total b — > ccs yield. The first number corresponds to x = 0.6 and the 
second to x — 1. Moreover, if our assumptions about the decay modes of D", D*', 
D*", D sl , D' sl , D* and D*' are justified, they lead predominantly to the "wrong" 
sign .D's. The yield of the "wrong" sign .D's then comes out to be 13 to 19%.^ The 
yield of the "right sign" .D's is estimated to be 10 to 12%; about 4% comes from 
the modes D^*'D^*\ and about the same from the modes with r\ in the final state. 
Thus, not only we confirm, exploiting a different line of reasoning, the conclusion 
of Ref. Q regarding the abundance of the "wrong" sign D's, but we encounter a 
menace of "overshooting": the yield of the "wrong" sign D's naturally tends to be 
too high. 

Under the circumstances it is tempting to calibrate theoretical predictions by 
using the CLEO and ALEPH data on the "wrong" sign D's, Br (B -> D^D^K) = 
10 ±4%. Combining this result with our estimates we see that the value of x ~ 0.6 is 
somewhat preferred. If so, the theoretically expected branching ratio of the "wrong" 
sign D's is about 13%, that of the "right sign" D's is about 10%, and the total 
branching ratio of the ccs channel is about 23%. Thus, we see that our model of 
saturation of inclusive data with two- and three-particle exclusive modes is perfectly 
consistent with the duality-based prediction of Ref. 0, ~ 23% or slightly higher. It 
is quite remarkable that the exclusive estimates presented above do not contradict 
the quark-hadron duality. This is not too surprising, of course, since in modelling the 
resonance yields we respected the constraints imposed by QCD plus factorization. 
The accuracy of this conclusion is not high, though, at the level of 30%. 

4 We assume here, rather arbitrarily, that the excited mesons from the be bracket are 

produced at the level of 2%, and this yield is split equally between the "wrong" sign D's and D s . 
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As an experimental confirmation of our approach the DD s r] + ... yield will have 
to be seen at the 4% level. 

The predictions for the absolute rates made above do not pretend to have too 
high an accuracy; suffice it to mention that the rates are proportional to /|> , 
and this quantity carries a noticeable uncertainty. However, it seems certain that, 
contrary to naive expectations, the modes with D s (the "right" sign .D's) constitute 
a smaller fraction of the b — > ccs channel. This seems to be a rather general feature 
not depending on details of the approximations we rely on. The spectral density 
emphasizes larger invariant masses of the cs system, which then predominantly 
generates decay chains with strangeness and charm separated. 

In our calculations we employed factorization, which must be violated at some 
level. However, if the cross-talking between two currents cs and h — > c appears 
to be strong, then one must expect dissolving D s into separate charm and strange 
particles from the very beginning. 

It is worth noting that the relative yield of is more sensitive to the strange 
quark mass, than the overall ccs yield. 

We point out that, although our calculations do not seemingly suggest a promi- 
nent mechanism for differentiating the hadronic widths of B and A& in the h — > ccs 
channel, where the calculations naively go in closest analogy to the one presented 
above, there still is a potential of locating the origin of differences. The standard 
arguments for factorization based on 1/N C counting rules are not necessarily appli- 
cable for the decays involving heavy baryons, in particular when the velocities of 
the final particles are small. This theoretical possibility must be carefully explored. 
On the other hand, a detailed experimental information about Ab hadronic width 
and, first of all, the value of n c , can provide us with an invaluable direct input. 
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